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Recently several lattice collaborations have studied the scale dependence of the coupling in the- 
ories with different gauge groups and fermion representations using the Schrodinger functional 
method. This has motivated us to look at the convergence of the perturbative step scaling to its 
continuum limit with gauge groups SU(2) and SU(3) with Wilson fermions in the fundamental, 
adjoint or sextet representations. We have found that while the improved Wilson action does re- 
move the linear terms from the step scaling, the convergence is extremely slow with the standard 
choices of the boundary conditions for the background field. We show that the situation can be 
improved by careful choice of the boundary fields. 
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1. Introduction and motivation 

Measuring the scale evolution of the coupling in gauge theories with matter fields in higher 
representation has been studied intensively over the past few years [j|]. It has been shown in [§,[3], ||] 
that for SU(2) and SU(3) with fundamental fermions the lattice step scaling function converges 
rapidly to its continuum limit when one uses &{a) improved Wilson-clover action. This can also 
easily be seen in the left panel of figure [jj where the fermionic contribution, normalized to the 
contunuum value, is shown. However, the situation changes when one considers fermions in the 
higher representations, which can be seen from the right panel of figure |]; see also [0]. 

Clearly the fermionic step scaling function for higher representations has large @{a 2 ) con- 
tributions, which are absent in the step scaling for the fundamental representation fermions. In 
the following we will study the effect of the Schrodinger functional boundary conditions on the 
convergence of step scaling [Q]. 



1.7 




Figure 1: Leading perturbative finite volume corrections to the continuum normalized fermion contribution 
of the lattice step scaling function with standard Schrodinger functional boundary conditions (described 
below). Fundamental fermions on the left panel and higher representation fermions on the right. 

2. Theoretical background 

We use the standard Wilson-clover lattice action 

S = S G + S F + 8S G , b + S gf + S FP , (2.1) 

with the usual single-plaquette action, So- The effect of higher representations appears in the 
Wilson fermion action, Sp, where the parallel transporter appearing in the gauge covariant lattice 
derivatives must be transformed into appropriate representation. 

With perturbative calculations, one must fix the gauge and this is reflected in the action via 
terms S g f and S^p. The specific form of these terms is only important while calculating the con- 
tribution of the gauge sector to observables. Since we will mainly be focused on the fermionic 
contribution, we refer to the original article ^ for details. 

For periodic boundaries, the 6(a) improvement is obtained by nonperturbatively tuning the 
coefficient c sw of the Sheikholeslami-Wohlert -term H^], 
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which we assume to be contained in the Wilson fermion action Sp- For the perturbative analysis 
described here we set c sw = 1. Also, for the fermion field "twist" in the spatial directions [|lC|]. 
\lf{x + Lk) = e l9k y(x), we use 6 = n/5 throughout. 

We wish to apply the Schrodinger functional scheme, i.e. introduce a constant background 
field by setting boundary conditions for the gauge fields at times 7 = and T = L, while retaining 
periodic boundaries in the spatial directions. This background field introduces further G(a) contri- 
butions which are cancelled by the boundary counterterms 8Scm in the action (|Q|). We will now 
describe the analysis of their effects in more detail. 

The boundary fields used for SU(2) are 





12 3 

; ; ; 



where 



<h = -n> <h = n, ti = n-p, ¥i = p-n- 

The standard choice for the angles are r\ = ? and p = n [11|. For SU(3) we used 
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(2.5) 



-T7(v + 2) + P> 
-0 2 + 2p. 



(2.6) 
(2.7) 



The standard choice for the angles are T] = 0, p = ? and v = [12]. 

With the boundary matrices Q and Cl from (2.3) or (2.5), we end up with a background field 
of the form 

B = 0, B k = (xtiC k + (L-xo)C k )/L, k = 1,2,3. (2.8) 
Using the path integral representation Z of the Schrodinger functional the effective action is 

r = -inz = g 2 r +r 1 + ^%2). (2 .9) 

The running coupling can be determined by studying how the system reacts to the change of the 
background field. Using the effective action, one defines the running coupling and obtains its 
perturbative expansion as 



8 2 (L) 



where 



dr /dri 
dF/dri 



Pi(L) 



d(r +ri)/d7] 

dTi/dn 



:gl +Pl (L)g 4 O +0(gl), 



(2.10) 



p h0 (L)+N F pi,i(L) 



dr /dr! 

and the quantity p\ has been split into gauge and fermionic parts, p\ o and p\\, respectively. 



(2.11) 
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The ff(a) errors arising from the fixed boundary conditions at times t = and t = T are 



removed by the countertems [13] 

1 



8S G j, = - I (c 1 -l)£Tr[l-[/( Av )] + -(Q-l)£Tr[l-[/(p f )], 
z <?o ps So P , 

8S F , h = a\c s - + &M +a\c t - - & t {x)]. 



(2.12) 
(2.13) 



For the precise form of the operators O sf and O s t we refer to [13]. Together with the bulk im- 
provement coefficient c sw , tuning the coefficients c s ,c t ,c s ,c t to their proper values we can remove 
all &(a) errors. 

For the electric background fields which we consider here, the terms proportional to c s do 
not contribute. Also, if we set the fermionic fields to zero on the boundaries, the counterterm 
proportional to c s vanishes. However the two terms proportional to c t and c t remain non- zero. The 
c t term corrects the mass of the fermions at times T = a and T = L — a and the c t term changes the 
weight of the time-like plaquettes on the boundary. In this work we set c t = 1 and c t to its one loop 
perturbative value [0, |], ^, [l2|]. This gives us one loop 6(a) improvement to the lattice step scaling 
function. 

The step scaling function and its perturbative expansion to one loop are 



Z(u,s,L/a) = g 2 (go,sL/a)\ g 2 {goL/a)=u , 

= M + [Ii ;0 (s,Va)+ £ i ! i( s ' L /«)NF]" 2 + 



(2.14) 
(2.15) 



Using the perturbative expansion of the coupling from fl2.10| ) in the perturbative formula for the 
step scaling function, we obtain for the one loop coefficients the following useful formulas: 



L L0 (s,L/a) = pi )0 (sL)-pifi(L), 
Ei,i (j,L/a) = pi,i(sL)-pi t i(L). 



We also introduce the variable 



Si 



I M (2,Z» Ei,/(2,L/a) 



i = 0,l, 



(2.16) 
(2.17) 



(2.18) 



(Ti,,- (2) 2fc 0li ln2 

which is the ratio of the perturbative step scaling and its continuum limit. This variable is useful in 
illustrating the convergence of the step scaling. In equation (2. 18) we used 

6o,o = 1 1AW(487T 2 ), 2>o,i = -7k/(127T 2 ), (2.19) 

which are the one loop coefficients of the perturbative beta function. 

The choice of the boundary fields in (14) and (2/7) is in no way unique. In fact one can 
choose the form of the boundary fields and the values of the angles T7 , p and v quite freely. The 
only limitation is that the fields and 0' have to belong to the so called fundamental domain. This 
consists of all the boundary fields that satisfy the equations 



01 < 02 < ... < (j>„, 



|0j — 0/| < 271, for all i,j, 



N 

I> = 0. 

i=i 



(2.20) 



4 



Convergence of step scaling 



Tuomas Karavirta 



Boundary fields of this type lead to a unique (up to a gauge transformation) minimal action, as has 
been shown in 



We take the fields and 0' that were introduced in equations (2.4) and (2/7) and choose 



variables T\ and p to be free parameters. The equation ( |2.20| ) then gives us the range of allowed 
values for rj and p. 



3. Numerical results 

Here we will present absolute errors of the the step scaling functions = 10, s = 2) as a 

function of v\ and p for the adjoint representations of SU(2) and SU(3) and the sextet representation 
of SU(3). The plotted variable thus is |5i i — 1|. 




Figure 2: Fermionic part of the Lattice step scaling function for SU(2) with adjoint fermions at L = 10. 
Optimal choice is p = % and f] = 5 . 

Consider first SU(2) with two adjoint Dirac fermions. There are two darker regions in the 
figure ^, indicating the areas where the discretization errors are the smallest. These two values of 
p, 7r/2 and 37r/2 are actually equivalent. The 77 dependence is weak, and v\ can be selected from 
the values within the fundamental domain quite freely (this turns out to be true also for SU(3)). 
However the value r\ = p/2 must be excluded, since the step scaling function diverges at that 
point. 




Figure 3: Fermionic part of the Lattice step scaling function for SU(3) with adjoint fermions at L = 10. 
Optimal choice is p = | and TJ = — f . 



5 



Convergence of step scaling 



Tuomas Karavirta 



Then turn to SU(3) with two Dirac fermions either in the adjoint or the sextet representation. 
Figure || shows also two darker areas for SU(3) with adjoint fermions. More detailed analysis 
shows, that in the p = n/6 area the &(a 2 ) effects are smaller. In figure [|, which shows the step 
scaling function for SU(3) sextet, there is only one dark region at p = 677T/150. We also found out 
that, these optimal values of v\ and p are independent of L. Change in L only modifies the scale of 
the errors. 




-0.8 -0.6 -0.4 -0.2 0.2 

T) / 71 



Figure 4: Fermionic part of the Lattice step scaling function for SU(3) with sextet fermions at L = 10. 
Optimal choice is p = and rj = — $. 




0.1 

a/L 
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Figure 5: Fermionic part of the Lattice step scaling function for higher representation fermions with the 
new boundary conditions. 



We also want to be sure that the value of the step scaling parameter s does not have an effect 
on the convergence of step scaling. While the convergence does improve slightly as s is increased, 
the change is still quite small and not enough to cure the higher order lattice effects present in the 
step scaling functions with the old boundary conditions 1 . We have also checked that s has no effect 
on the preferred values of r\ and p . 

One has to keep in mind that changing the boundary conditions could in principle also affect 
the convergence of the gauge part of the step scaling function. This happens to be true for SU(2), 



'From now on we will use the terms "old" and "new" boundary conditions. The old boundary conditions refer to 
the conditions presented in equations ( J2.4| ) and ( 2/7 ) with the standard choice of the parameters p and r\ . With the term 
new boundary conditions we refer to the optimal choice of the parameters p and r\ for a specific symmetry group and 
fermion representation as stated in this section. 
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but fortunately the effects are roughly a hundred times smaller. Thus one can freely choose the 
boundary conditions that are optimal for the fermionic part without compromising the convergence 
of the step scaling function. In fact the convergence is even faster if one uses the new boundary 
conditions that are optimized for the adjoint fermions. We expect similar behavior from the gauge 
sector of the other gauge groups. 

4. Conclusions 

We have provided a cure to the slow convergence of the lattice step scaling function with 
higher representation fermions. Our method has the drawback, that in simulations with higher 
representation fermions and new boundary condition the signal for measuring the scale evolution 
of the coupling constant can be smaller than it is with the old boundary conditions. However one 
has to get rid of the higher order lattice artifacts, to get reliable results. Other solutions to this 



problem exist for SU(3) [14p. A combination of these methods might provide a way to get the 



smallest possible lattice artifacts while retaining also the signal optimally high. 
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